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The Symanzik analysis of correlators in lattice QCD with maximally twisted Wilson fermions 
reveals that there exist cutoff artifacts which tend to become large as the quark mass gets small. 
We show that these effects can be reduced to a negligible level by either introducing the clover 
term in the action or, as already suggested in the literature, by a suitable choice of the critical 
mass. Recent simulation data support these conclusions. 
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1. Introduction and main results 

Although in lattice QCD at maximal twist (Mtm-LQCD) 0(a) discretization effects (actually 
all O(a 2k+l ),k>0 effects) are absent or easily eliminated |3L 0], it turns out that correlators 
are affected by dangerous artifacts of relative order a 2k , k > 1, which are enhanced by inverse 
powers of the (squared) pion mass, as the latter becomes small. In fact, when analyzed in terms of 
the Symanzik expansion, lattice expectation values exhibit, as m 2 K —> 0, what we will call "infrared 
(IR) divergent" cutoff effects with a behaviour of the form 

L cont r / q^ 1 \ - 

(O) ={0) 1+0 , 2k > h > 1 (k,h integers), (1.1) 

where we have assumed that the lattice correlator admits a non-trivial continuum limit. Powers of 
Aqcd required to match physical dimensions are often understood in the following. 

We shall see that artifacts of the type ( |1.1| ) are reduced to terms that are at worst of order 
a 2 (a 2 1 'm^) fc_1 , k > 1, if the action is 0(a) improved a la Symanzik, or alternatively the critical 
mass is chosen in some "optimal" way. 

The idea that a suitable definition of critical mass exists which can lead to a smoothing out 
of chirally enhanced lattice artifacts or perhaps be of help in getting improvement was already put 
forward in the context of chiral perturbation theory in refs. ^ and [g], respectively. 

An important consequence of our analysis is that the strong (order of magnitude) inequal- 
ity m q > oAqcpj, invoked in ref. ^ can be relaxed to the weaker relation m q > a 2 Ag CD , before 
large cutoff effects are possibly met while lowering the quark mass at fixed a. The works of 
refs. [Q, and most recently refs. [|[], which are based on lattice chiral perturbation theory, 
lead to essentially equivalent conclusions about cutoff effects in pion quantities in the parameter 
region m q > <3 2 Aq CD . They also yield interesting predictions on the possible Wilson fermions phase 
scenarios []9[ [T(|] and results, when m q is of order a 2 or smaller. 

A thorough discussion on the effectiveness of Mtm-LQCD in killing 0(a) discretization errors 
and the ability of the optimal choice of the critical mass in diminishing the magnitude of lattice 



artifacts at small quark mass can be found in [ |11| , |12| ] and in the work of refs. [|13j, 140. As for Mtm 



LQCD with clover-improved quark action, the promising quenched tests presented some years ago 
in [15] have been recently extended in Jl6| ] down to pion masses of 300 MeV or lower, confirming 



the absence of large cutoff effects. 

The outline of this presentation is as follows. In Section || we analyze the form of the Symanzik 
expansion of lattice correlators beyond 0(a) and explain why and how TR divergent" cutoff effects 
arise in this context. In section || we discuss two ways of killing all the leading "IR divergent" cutoff 
effects and we describe the structure of the left-over "IR divergent" terms. Finally in Section |] we 
collect some remarks on the peculiar structure of the lattice artifacts affecting lattice hadronic 
energies and in particular pion masses. Conclusions can be found in Section [|. 

2. Symanzik analysis of "IR divergent" cutoff artifacts 

The study of cutoff artifacts affecting lattice correlators in Mtm-LQCD can be elegantly made 
in the language of the Symanzik expansion. A full analysis of cutoff effects beyond 0(a) is of 
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course extremely complicated. Fortunately it is not necessary, if we limit the discussion to the 
terms that are enhanced as the quark mass m q is decreased. 

• The Symanzik LEEA of Mtm-LQCD - The expression of the fermionic action of Mtm-LQCD 
in the physical quark basis is given in |3||. The low energy effective action (LEEA), Ss ym , of the 
theory can be conveniently written in the form 



/ CO CO 
d 4 y Jgf 4 (y) + £ a 2k+l £ 4+2k+l (y) + £ a 2k £ 4+2k (y) 
k=0 k=l 



(2.1) 



where S£ 4 = ^jtr(F -F) + y{y D + m q ) \ff is the target continuum QCD Lagrangian. Based on 
the symmetries of Mtm-LQCD a number of interesting properties enjoyed by 5s ym can be proved 
which are summarized below. 

1. Lagrangian densities of even dimension, £ 2k , in eq. ( |Q| ) are parity-even, while terms of odd 
dimension, £ 2k+ \, are parity-odd and twisted in iso-spin space. Thus the latter have the quantum 
numbers of the neutral pion. 



2. The term of order a in eq. ( |2.6[ ), £5, is given by the linear combination 



^5 = S 5:SW £ 5:SW + S 5 m 2 £ 5 m 2 + 8$ ie £ 5:e , (2.2) 

is,sw = ■ F}iy 5 i: 3 Y, 4m? = ™ 2 q W, he = Aqcd Y'WsV, ( 2 - 3 ) 

where the coefficients ds^w, 8 5m i and 85^ are dimensionless quantities, odd in r. The operator £5^ 
arises from the need to describe order a uncertainties entering any non-perturbative determination 
of the critical mass and goes together with £$ t sw- Both £5^ and £5^ could be made to disappear 
from (^jj) by introducing in the Mtm-LQCD action the SW (clover)-term [|l7|] with the appropriate 
non-perturbatively determined csw coefficient [|lj|] and at the same time setting the critical mass to 
its correspondingly 0(a) improved value. 

3. Higher order ambiguities (k > 1) in the critical mass, which will all contribute to Jz? dd> are 
described by terms proportional to odd powers of a of the kind 

a 2k+l 8 4+2k+he £ 4+2 k+i,e = a 2k+l 8 4+2k+l , e (A QCD ) 2k+2 yiy 5 ^y . (2.4) 

• Describing Mtm-LQCD correlators beyond 0(a) - We are interested in the Symanzik descrip- 
tion of the lattice artifacts affecting connected expectation values of «-point, multi-local, multiplica- 
tive renormalizable (m.r.) and gauge-invariant operators 0(x\,x 2 ,. . . ,x n ) = IX/=i Oj{ x j) = 0(x), 
ii / ... ^4, which we take to have continuum vacuum quantum numbers, so as to yield 
a non trivially vanishing result as a — > 0. In order to ensure automatic 0(a) improvement [Q] we 
shall assume that O is parity invariant in which case its Symanzik expansion will contain only even 
powers of a. Schematically we write 

(0(x)) L = (\0(x) + A odd <9(x) +A even O(i)] e -^ 4 >'^M +y -(>')l) C ° nt , (2.5) 

m q m, 

CO CO 

-^odd = £ a 2k+l £ 4+2k+ i , -Steven = £ ^ ^4+2* • (2-6) 

k=0 k=\ 

The operators A Q ddO (A even O) have an expansion in odd (even) powers of a. They can be viewed 
as the «-point operators necessary for the on-shell improvement of O [jl9[ [T^] . 
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• Pion poles and "IR divergent" cutoff effects - Although a complete analysis of all the "IR 
divergent" cutoff effects is very complicated, the structure of the leading ones (h = 2k in eq. (PI)) 
is rather simple, as they only come from continuum correlators where 2k factors / d 4 y Jz? dd(.y) are 
inserted. More precisely the leading "IR divergent" cutoff effects are identified on the basis of the 
following result [|]]. 

In the Symanzik expansion of (0(x))\^ n at order a 2k (k > 1) there appear terms with a 2k- 
fold pion pole and residues proportional to |(£2|«£f o dd|ft (0))| 2 *, where (H| and |?r (0)) denote 
the vacuum and the one-7T° state at zero three-momentum, respectively. Putting different factors 
together, each one of these terms can be seen to be schematically of the form (recall Jz? dd = 0(a)) 

cont 



] \2k -i cont 

-=-) (^) 2k ^[O;{7i°(0)} 2 k} , (Hl^oddlAO)) 



77? ^ 
rn 7t 



(2.7) 



where we have generically denoted by .J{\0\ {^(O)^] the 2fc-particle matrix elements of 0, with 
each of the 2k particles being a zero three-momentum neutral pion. 

Less "IR divergent" cutoff effects (those with h strictly smaller than 2k in eq. (|1.1|)) come either 
from terms with some extra / d 4 yJ£' even (y) insertions or from contributions of more complicated 
intermediate states other than straight zero three-momentum pions or from both. In the first case 
one gets extra a powers (not all "compensated" by corresponding pion poles), while in the second 
one loses some 1 / m\ factor. 



It is important to remark that the appearance of pion poles like the ones in eq. fl2.7p in no way 
means that the lattice correlators diverge as m q — > 0, but only that the Symanzik expansion we have 
employed appears to have a finite radius of convergence (on this point see the remarks of ref. [^]). 

3. Reducing "IR divergent" cutoff artifacts 

Recalling that Jz? dd = a £5 +0(a 3 ), the previous analysis shows that at leading order in a the 
residue of the most severe multiple pion poles is proportional to | (Sl\£s |tt°(0)) \ 2k . It is an immediate 
conclusion then that the leading "IR divergent" cutoff effects can all be eliminated from lattice data 



if we can either reduce £5 to only £ 5 m i in (2.2) or set ^ n to zero 



Improving the Mtm-LQCD action by the SW-term - The obvious, field-theoretical way to 



eliminate £5 from the LEEA of Mtm-LQCD consists in making use of the 0(a) improved action [ 17 
HHH- In this case lattice correlation functions will admit a Symanzik description in terms of a 
LEEA where the operators £s,sw and £5^ are absent, and £5 is simply given by £ 5 m 2. The left-over 
contributions arising from the insertions of £ 5m 2 in (0}|™ nt yield terms that are at most of order 
(am 2 /m 2 K ) 2k ~ (am q ) 2k , hence negligible in the chiral limit. It is instead the next odd operator in 
the Symanzik expansion, a 3 £7, which comes into play. 

A detailed combinatoric analysis based on the structure of the non-leading "IR divergent" 
cutoff effects [Q] reveals that the worst lattice artifacts left behind in correlators after the "clover 
cure" are of the kind a 2 (a 2 /m^) , k > 1. 

• Optimal choice of the critical mass - The alternative strategy to kill the leading "IR divergent" 
cutoff effects consists in leaving the Mtm-LQCD action unimproved, but fixing the critical mass 
through the condition 



lim ^K)= lim <n|j§f odd |7T°(0)} 

n„— >0+ m„— >0+ 



cont 

= 0. (3.1) 
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The meaning of (3.1) is simple. It amounts to fix, for k > 0, the order a 2k+{ contribution in the 



counter-term, M cr \f/ L iy5r3\j/ L , so that its vacuum to one-7T°(0) matrix element compensates, in the 
limit m q — > 0, the similar matrix element of the sum of all the other operators making up ^4+2/t+i - 

A concrete procedure designed to implement condition ( |3.1| ) in actual simulations was dis- 
cussed in ref. [Qj. It consists in determining the critical mass by requiring the lattice correlator 
a 3 J2 X (Vq(x)P 1 (0)) |£ (jco 7^ 0) to vanish in the chiral limit, where V 2 = y/^yl/Ais the vector cur- 
rent with iso-spin index 2 and P 1 = yjs^y the pseudo-scalar density with iso-spin index 1. In the 
continuum this correlator is zero by parity for any value of m q . On the lattice the breaking of parity 
(and iso-spin) due to the twisting of the Wilson term makes it non-vanishing by pure discretization 
artifacts, which have the form of a power series expansion in ^/m 2 . 

The important conclusion of the analysis presented in is that it is not necessary (nor possi- 
ble) to really go to m q — > 0. It is enough to have the critical mass determined by the vanishing of the 
above correlator at the current simulation quark mass, provided we stay in the region m q > a 2 . In 
these conditions we will have E,%{m q ) = 0(am 2 ) with all the leading "IR divergent" cutoff effects 
reduced to finite 0{a 2k ) terms. As for the subleading ones, a non-trivial diagrammatic analysis 
shows that the worst of them, left behind after the "optimal critical mass cure", are reduced to only 
a 2 (a 2 / m 2 ,)* -1 , k > 1, effects, just like in the case where the clover term is employed. 

4. Artifacts on hadronic energies and pion masses 

In the language of the Symanzik expansion discretization artifacts on hadronic energies are 
described by a set of diagrams where at least one among the inserted / Jz? dd factors gets neces- 
sarily absorbed in a multi-particle matrix element, with the consequence that it is not available for 
producing a pion pole. As a consequence, at fixed order in a, the most "IR divergent" lattice cor- 
rections to continuum hadronic energies contain one overall factor l/m 2 less than the leading "IR 
divergent" cutoff effects generically affecting correlators. For instance, to order a 2 the difference 
between lattice and continuum energy of the hadron a„ reads [Q] 



A£« n (q) 



^(^|£ 5 |7r (0))(7r°(0)a ;1 (q)|^|a„(q))\ , n , 2 ^ 



-Re ( > v " v ^ 1 +0{a 



(4.1) 



where 0(a 2 ) denotes "IR finite" corrections. It should be noted that this "IR divergent" lattice 
artifact is reduced to an "IR finite" correction after anyone of the two "cures" described in Sect. ||. 



Specializing the formula ( |4.1| ) to the case of pions, one obtains the interesting result that the 
difference between charged and neutral pion (square) masses is a finite 0(a 2 ) quantity even if the 
critical mass has not been set to its optimal value or the clover term has not been introduced. The 



reason is that the leading "IR divergent" contributions shown in ( |4. 1[ ) are equal for all pions (as 



one can prove by standard soft pion theorems [20]), hence cancel in the (square) mass difference. 



This conclusion is in agreement with detailed results from chiral perturbation theory (see refs. JTo| ] 
and []|]), as well as with the first numerical estimates of the pion mass splitting in Mtm-LQCD 

5. Conclusions 

When analyzed in terms of the Symanzik expansion, lattice correlators in Mtm-LQCD show 
"IR divergent" cutoff effects which tend to become large as the quark mass gets small. Extending 
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the works of refs. [Q, ^J], we have shown that, not only if the critical mass is chosen in some 
"optimal" way but also if the action is clover improved, such lattice artifacts are strongly reduced to 
terms that are at worst of the type a 2 (a 2 /m|) , k > 1. The latter result implies that the continuum 
extrapolation of lattice data is smooth at least down to values of the quark mass satisfying the order 
of magnitude inequality m q > a 2 Ag CD . 
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